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The approximations are based on spectral properties of the noise. They allow one to consider the 

noise as the result of fractional integration/differentiation of a white Gaussian noise. We study 

correlation properties of the approximation to fractional Gaussian noise and point to the 

peculiarities of persistent and anti-persistent behaviors. We also investigate self-similar properties 



O . TT 

Q , of the approximation to fractional Brownian motion, namely, "r laws" for the structure 



function and the range. We conclude that the models proposed serve as a convenient tool for the 
natural processes modelling and testing and improvement of the methods aimed at analysis and 

interpretation of experimental data. 
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I. INTRODUCTION. 

Fractional Brownian motion (fBm) is a random process with stationary self-similar in- 
crements, which are normally distributed and have an infinite span of interdependence. The 
general theory of random processes with stationary self-similar increments was developed 
by Kolmogorov ||T| as far back as 1940. Fractional Brownian motions were introduced by 
Mandelbrot and van Ness as a (relatively) simple family of random functions "that could 
in some way be expected to be "typical" of what happens in the absence of asymptotic 
independence " 0] . In this paper the authors also have introduced fractional Gaussian noise, 
that is, the stationary process, which is the derivative of the "smoothed" fBm. Three classes 
of examples have moved these authors to study fractional motions, namely, (i) various eco- 
nomic time series, (ii) 1// , or flicker, noises in electronics and solid state physics, and (iii) 
the tasks of hydrology, see references in 0. Over past 30 years the usefulness of fBm's was 
conflrmed repeatedly, and the number of phenomena, which can be viewed as examples of 
fBm's grows. There are many papers on this subject; here we only mention the monograph 
0], which contains a lot of natural examples, the review 0], which deals with the problems 
of biology and physiology, and the recent monograph |^ , in which statistical properties of 
economical and flnancial time series are discussed. We also that the description of fractional 
noises in terms of probability densities requires the use of fractional calculus, which has an 
ever growing interest among physicists and economists [^. 

The ubiquity of fBm's has posed the task of developing approximations, which allow 
one to simulate fBm's and study their properties in order to develop effective methods for 
experimental data processing and for the purposes of numerical modelling. Soon after the 
flrst paper on fBm Mandelbrot and Wallis have carried out simulation experiments con- 
cerning fractional Gaussian noise . They have proposed two approximations to the noise 
called "Types 1 and 2 approximations". Later on, Voss have proposed an empirical algo- 
rithm PI suited for construction of not only one-dimensional fractional Brownian functions 
but also fractional Brownian surfaces and volumes. However, we believe, that the ways of 



constructing fractional noises and motions are not exhausted, and various simulation models 
are needed, each of them may appear to be useful when studying some particular problem. 
In this paper we propose a simple model of fBm based on exploitation of spectral prop- 
erties of fractional Gaussian noise. At first we remind some features, which are essential for 
our simulation and studies. 

II. SOME PROPERTIES OF FRACTIONAL GAUSSIAN NOISE AND 
FRACTIONAL BROWNIAN MOTION. 

Below we use the same notations as in original papers by Mandelbrot and van Ness 
and by Mandelbrot and Wallis [0. An ordinary Brownian motion with zero mean and unit 
variance is denoted by B{t), and fBm by Bij{t). The latter is constructed as a moving 
average of dB(t), in which past increments of B{t) are weighted by the kernel (t — t')^^^^"^ , 
where if is a parameter satisfying < if < 1. For the ordinary Brownian motion H = 1/2. 
For H lying between 0.5 and 1 the process is called persistent, whereas for H lying between 
and 0.5 it is called anti-persistent. The introduction of fBm is motivated by the self-similar 
properties of its increments, which lead to the "r-'^ laws" for the structure function and for 
the range, that is, 

sTir) = ({Bnit + r) - B^it))')'^' = r^vT (1) 



(an explicit form of Vu is unessential here), and 

R{t,T)= sup \Bn{s)-Bu{t)\- 

t<s<t+r 

- inf [BH{s)-BH{t)]^ 

t<S<t+T 

= i?(0,r) = r^i?(0,l). (2) 

The properties of self-similarity of approximation to fBm will be checked below. 
It is known that the sample paths of fBm are nowhere different iable. Such irregular 
functions are best characterized by the Holder exponent or by recently introduced local 



fractional derivative |^. For the purposes of our paper we note that in apphcations one 

always deals with the "smoothed" random process, Buit^d) , which may be defined as 

t+s 
BHit,6) = - J BH{s)ds , (3) 

t 

where S is the smallest time interval, which is physically relevant for the problem of interest. 
Buit.d) has a derivative, 

BH{t,5) = -^[BH{t + 5)-BH{t)], (4) 

which is called fractional Gaussian noise. It is a stationary process with the covariance, 
which may be easily determined, 

C{T,5) = {B'^{t + r,5)B'^{t,5)) 

If r ^ (5 , it follows from Eq.(5) that 

C{T,6)(x{2H-iy-\ (6) 

and thus, the covariance has the same sign as if — 1/2 and has a very different behavior in 
a persistent case and in anti-persistent one. The analysis of Eqs.(5), (6) shows that for the 
former case C(r, 5) is positive for all r, and 

oo 

/ Ch{s., 5)ds = cx), 



whereas for the latter C(r, 6) changes sign once from positive to negative at a value of r 
proportional to 6, and one has 



f CH{s,6)ds = 0. 





The next property, which serves as the basic one at the initial stage of our simulation, is a 
power-law behavior of the spectrum of fractional Gaussian noise. 



G{f; 6)=AJ Ch{s, 6) cos{27t fs)ds oc /^"^^ (7) 



at f6 « 1. 

The spectrum (7) is also quite different for persistent and anti-persistent cases. For the 
former the main part of the noise energy is concentrated in a low -frequency region, whereas 
for the latter the spectral intensity grows towards high frequencies. 

In numerical simulation one deals with a discrete-time fractional noise, when t is a 
positive integer and the noise is defined by the sequence of increments Bnit + 1) — Buit) 
. The correlation function is defined by Eq.(5) with 5=1. The Fourier transform of the 
discrete-time noise is defined between plus and minus the Nyquist critical frequency /c = 0.5. 

The power-law dependence of the spectrum at small / allows one to simulate fractional 
Gaussian noise by fractional integration/differentiation of a white noise, with a subsequent 
construction of fBm. 

III. THE MODELS. 

We remind the relation between the Fourier transforms of the function X{t) and of its 
fractional integral/derivative Xy{t) of the v order: 

L if) ^ ^, (8) 

where z/ is positive in case of fractional integration and negative in case of fractional differ- 
entiation (to be more accurate, we say about left-side Liouville fractional integral/derivative 
on the infinite axis ||lO|). If X(t) is a white Gaussian noise, then the spectral density of 



Xi,{t) is, at least, formally. 



1 



G.^oc— , (9) 



which behaves as the spectrum of fractional Gaussian noise at small /, see Eq.(7), if we set 
V = H — 1/2, —1/2 < V < 1/2 . Thus, we use the following way for simulation: 



(i) taking white Gaussian noise X{t), t is an integer, we multiply its Fourier transform 
by /-^ -1/2 <u< 1/2. 

(ii) making an inverse Fourier transform, we get X^{t), which is supposed to approximate 
fractional Gaussian noise with the index H = +1/2. 

(iii) The process X*(t) , which is supposed to approximate fBm, is defined by 

x:(t) = ^x.(t). (10) 

r=l 

IV. NUMERICAL RESULTS. 

The results of numerical simulation and analysis are shown in Figs. 1-5. 

At the left of Fig.l we show typical examples of the samples of the process X,y{t) for 
(a) u = 0.4, (b) u = 0(white Gaussian noise), and (c) z/ = —0.4. The stationarity of the 
samples was checked by the methods of series and inversions |TT[. The figures demonstrate 



clearly the prevalence of low-frequency components for the case (a), and the prevalence of 
high-frequency components for the case (c). At the right of Fig.l the solid lines depict 
normalized theoretical correlation functions, see Eq.(5), with 6 = 1, for (a) H = 0.9, (b) 
H = 0.5, and (c) H = 0.1, whereas the normalized correlation functions 

{XUt)XUt + r)) 

are shown by black points for the same values of as at the left of the figure. The results 
demonstrate good agreement between correlation function of simulated process with the pa- 
rameter z/ and theoretical correlation function for the fractional noise with the corresponding 
index H = u + 0.5. We also see clear difference between correlation functions for the persis- 
tent and anti-persistent cases. This difference is the manifestation of the property, which at a 
qualitative level can be formulated as follows: in the persistent random process the available 
tendency is supported, whereas in the anti-persistent process the opposite tendency prevails 
M M. We also note that the correlation functions indicated by black points are estimated 
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from the samples shown at the left with the use of time averaging. The complementary 
averaging over a set of realizations diminishes scattering of the points. 

It is known, that the wavelet transform is well-suited for analyzing structure of non- 



stationary processes |T2[ . In particular, it allows one to study the behavior of processes at 
different scales and hence it is sometimes called as a mathematical microscope. The wavelet 
transform T{a,t) of the process X*{t) is written as 

1 r „Jt~t' 



T{a,t) = — fw{- ]x:{t')dt', 

\ a J \ a I 
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where W{t) is called analyzing wavelet. In Fig. 2 we demonstrate the wavelet transforms 
T{a,t) of the samples of X*{t) for (a) u = 0.4, and (b) u = —0.4. We use the "Mexican 
hat" as analyzing wavelet, 

W{t) = (l-t2)exp(-tV2). 

The low part of the wavelet corresponds to small a values (that is, to high frequencies), 
whereas the upper part corresponds to large a values (low frequencies). Below each of the 
wavelet the sample path is shown, which is subjected to wavelet transform. It is seen that 
the wavelet for anti-persistent case has a clearly visualized structure in a high-frequency 
region. On the contrary, the wavelet for persistent case demonstrate almost homogeneous 
background in a high-frequency region. The lines of local maxima demonstrate hierarchial 
structure, which reflects self-similarity of fBm. 

In this paper we do not intend to use wavelets for extracting quantitative characteristics, 
so. Fig. 2 has an illustrative purpose only. However, we note that the wavelet transform 
being applied to some known functions allows one to establish their multifractality |jl3|. 



We study numerically r-dependence of the structure function of the generated process 

xm , 

S{r) = {{X:{t + T)~X:{t)f)^r'^ (12) 

Since for fBm the "r^ law" is fulfilled, see Eq.(l), we expect that s is close to i/ + 0.5. In 
Fig. 3 s depending on u is depicted by points for the persistent processes and by crosses for 



the anti-persistent ones. The expected relation is indicated by dotted hne. We see that the 
numerical relation is well fitted by the expected line, however, some discrepancies appear 
for the strongly persistent and strongly anti-persistent processes as well. At the insets the 
structure function versus time delay r is indicated by black points in log-log scale. The 
slopes of fitted lines give the values of s for z/ = —0.4 (top inset) and v = 0.4 (bottom inset), 
respectively. 

We now turn to the study of the range of X*{t) . In the empirical rescaled range analysis 
, that is, at experimental data processing or in numerical simulation the range of the random 
process is divided by the standard deviation of its increments after subtraction of a linear 



trend, see Ref. [|T4|. This procedure, in particular, smooths the variations of the range on 
different segments of time series. As the result of the empirical rescaled range analysis of 
experimental data one gets the Hurst exponent of the process. For fBm the Hurst exponent 
thus obtained must coincide with the parameter H, see Eq.(2). Fig.4 demonstrates the 
application of the rescaled range analysis to the sample paths of X*{t) . The index u is 0.4, 
thus the value of the Hurst exponent, which we expect to get according our way of simulating 
fBm is 0.9. In Fig.4a the fluctuations of the range R (thin curve) and those of the standard 
deviation a of the increments (thick curve) are shown for the case when the total length 
of the sample is divided into 64 segments, each of =16 lengthwise. Below the variations of 
the ratio R/a are depicted. It is shown that fluctuations of the ratio are smaller than those 
of the range. This circumstance justifies the use of the ratio in the empirical analysis. In 
Fig. 4b the rescaled range versus time interval r is depicted by black points in log-log scale. 
The slope of the fitted line gives the Hurst exponent 0.83. 

In Fig. 5 the values of the Hurst exponent H,y obtained by applying the rescaled range 
analysis to the sample paths of X*{t) are shown by black points for the persistent processes 
and by crosses for the anti-persistent ones. The dots indicate the line H,^ = u + 0.5, that is, 
the relation, which we expect to get according our way of simulating fBm. It is shown that 
(i) fairly good agreement is for intermediate z/, (ii) deviations between the results for X*{t) 
and the theory for fBm grow to the ends of the interval of z/, and (iii) these deviations are 
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larger for the anti-persistent case than for the persistent one. 

V. DISCUSSION. 

Figures 3 and 5 allow one to conclude that in the persistent case there is a better 
agreement with the expected lines than in the anti-persistent case (with the length of the 
sample and discretization step being equal for both cases). In other words, in our model 
the low frequency errors are less important than the high-frequency ones. This conclusion 
is analogous to that concerning "Type 2 approximation" constructed by Mandelbrot and 
Wallis . This is not surprising since the way of its constructing in real space is similar to 
our way of constructing approximation in Fourier space. Moreover, we compare the results 
of our numerical simulations with those performed with the "Type 2 approximation" and the 
Voss algorithm ^ and found unessential differences between the results for the exponent of 
the structure function, see Fig.3, and for the Hurst index, see Fig. 5, for the three methods. 
We may conclude that our method is suited for performing simulation of fBm as well as the 
"Type 2 approximation" and the Voss algorithm are. In all three models the persistent case 
is simulated with higher accuracy than the anti-persistent one. Fortunately, the persistent 
case is prevalent in nature P], therefore, the method proposed in our paper can be used in 
various applications. At last, we believe that our method will be useful for direct constructing 
approximations to non-Gaussian fractional noises, in particular, fractional Levy noises. This 
topic will be the subject of forthcoming paper. 
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FIGURES 
FIG. 1. At the left: typical samples of the simulated noise X^it) for (a) u = 0.4; (b) v = 

(white Gaussian noise), and (c) i' = —0.4. At the right: normalized correlation functions of 

the simulated noises with the same u (black points). Solid lines indicate normalized theoretical 

correlation functions, see Eq.(5) with 6 = 1, for (a) H = 0.9, (b) H = 0.5, and (c) H = 0.1. 

FIG. 2. Wavelet transforms of the simulated process X* (t) for (a) z^ = 0.4, and (b) i^ = —0.4. 
Below each of the wavelets the sample paths subjected to wavelet transform are shown. 

FIG. 3. The exponent s in the structure function of the simulated process, see Eq.(ll), versus 
the index i^. Black points indicate persistent case, whereas the crosses indicate the anti-persistent 
one. The line s = u + 1/2 is depicted by dots. At the insets: the square root of the structure 
function versus time in log-log scale (black squares) for i/ = —0.4 (top inset) and i^ = 0.4 (bottom 
inset). The slopes of fitted lines give the values of the exponent s, see the main figure. 

FIG. 4. (a) Variations of the range of X* (t) (thin curve), of the standard deviation of Xi,{t) 
(thick curve), and of their ratio (below) at the different time segments, (b) Rescaled range versus 
time interval in log- log scale (black points). The slope of the fitted line gives the Hurst index 0.83. 

FIG. 5. Plots of the Hurst exponent Hiy versus i/ for the persistent case (black squares) and 
anti-persistent one (crosses). The line H^, = v + 1/2 is indicated by dots. 
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